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Abstract 

We present a revised interpretation of recent analysis of supernovae data. We evaluate the effect of the priors on 
the extraction of the dark energy equation of state. We find that the conclusions depend strongly on the Qj^ prior 
value and on its uncertainty, and show that a biased fitting procedure applied on non concordant simulated data 
can converge to the "concordance model". Relaxing the prior on VLm points to other sets of solutions, which are not 
excluded by observational data. 
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The existence and nature of dark energy is one of the 
most challenging issues of physics today. The publica- 
tion of high redshift supernovae discovered by the Hub- 
ble Space Telescope, bythe SCP collaboration and 
recently by Riess et al. Q , has been interpreted as agree- 
ment of the data with the so named AC DM "concor- 
dance model" (fi M ~ 0.3, Q A w 0.7, w=p/p = -1). We 
have reconsidered some conclusions in the light of our 
previous analysis of simulated data Q. 

Riess et al.|2| have selected 157 well measured SNIa, 
which they call the "gold" sample, a set of data we will 
use throughout this paper. Assuming a flat Universe 
(fix = 1) they conclude that: i) Using the strong 
prior Qm — 0.27 ± 0.04, a fit to a static dark energy 
equation of state yields -1.46< w <-0.78 (95%CL); 
ii) Looking at a possible redshift dependence of w(z) 
(using w(z) = wq + w\z), the data with the strong 
prior indicate that the region w\ < and especially the 
quadrant (wo > — 1 and w\ < 0) are the least favoured. 
They reject large time variation and are compatible with 
the concordance model. 

We have shown in that it is unavoidable to get 
some ambiguities when trying to fit a particular fiducial 
cosmology with a "wrong" model. This "bias problem" 
has been mentioned several times in the literature, see 
e '/ [3 00 In this letter, we explore the effect of the 
J]m prior on the determination of w(z). s = 

Following j3|, we assume a flat universe and keep the 
same paramctrisation of w[z) as in [2j, for the sake of 
comparison. We call 3-fit (4-fit) the fitting procedure 
which involves the 3 (4) parameters Mg, Qm and wq 
(Ms, £Im, and u>i), Ms being a normalisation pa- 
rameter (see y| for definitions and formulae). We have 
performed 3-fits and 4-fits and compared the results in 
different cases, varying the central value and the uncer- 
tainty on the Qm prior. 

TABLE I: Fit results obtained using the gold data from Q for 
various fitting procedures. The x 2 ls very stable, it is around 173 
(for 157 SNIa) for all procedures except for the 3-fit with the strong 
prior Q M = 0.27 ± 0.04 where \ 2 ~ 176. 



Fit 


fiju prior 




w 




3-fit 


no 


0.48 ± 0.06 


-2.2 ±0.95 


/ 


3-fit 


0.27 ±0.2 


0.45 ± 0.07 


-1.9 ±0.73 


/ 


3-fit 


0.50 ±0.2 


0.48 ± 0.06 


-2.3 ±0.94 


/ 


3-fit 


0.27 ± 0.04 


0.28 ± 0.04 


-1.0 ±0.15 


/ 


3-fit 


0.50 ± 0.04 


0.49 ± 0.03 


-2.5 ±0.77 


/ 


4-fit 


no 


0.48 ± 0.20 


-2.2 ± 1.34 


0.12 ±23 


4-fit 


0.27 ±0.2 


0.35 ±0.18 


-1.6 ±0.80 


1.74 ± 1.3 


4-fit 


0.50 ±0.2 


0.49 ± 0.20 


-2.6 ± 1.20 


1.60 ±18 


4-fit 


0.27 ± 0.04 


0.28 ± 0.04 


-1.3 ±0.26 


1.50 ±0.84 


4-fit 


0.5 ± 0.04 


0.49 ± 0.04 


-2.6 ± 1.40 


0.95 ± 10 



Applying no prior or the strong prior on £Im (lines 1, 
4 and 9 of the Table), we recover the results obtained by 
Riess et al.[2j. Nevertheless, some interesting points can 
be underlined: 



• With no prior or a weak prior on £Im, the preferred 
flju values are always greater than 0.3. 

• Without any assumption on U,m nor w\, the error on 
Q.M is close to 0.2 (line 6 of Table I). 

• Changing the central value of the f2j\/ prior leads to a 
change in the wq values of more than la. The u>o values 
are strongly correlated to J1m and are thus always smaller 
than the ACDM value, when the strong prior on Qm is 
relaxed. \ 2 1S very stable but the correlation matrix can 
vary a lot for the 4-fits and the (wo,wi) solution. 

• If the Qm prior is strong, the conclusion on wq depends 
on the prior value : for Om=0.27, Wo is forced to values 
compatible with -1, in particular for the 3-fit and the 
errors are strongly reduced. For S1m=0.5, wq is more 
negative and the errors are significantly larger. 

• The only cases where "reasonable" errors can be found 
on wi occur for £Im around 0.3. 

To illustrate these points, Figurc^Jshows the results in 
the (£Im,wq) plane for the 3-fits(left) and the 4-fits(right), 
using no prior on fijvf or two strong priors with the two 
central values: 0.27 and 0.5. As expected the contours 
strongly depend on the procedure used to analyse the 
data. For instance, the 95%CL contours for the two 
strong prior cases are disconnected. However, we note 
that f^M < 0.6 is valid for all procedures, hence it is one 
of the strong conclusions from present SN data. 
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FIG. 1: 95%CL contours for 3-fits(left) and 4-fits(right) with no 
prior on fijvf (plain) and two strong priors flu = 0.27(0.5) ± 0.04 
(dashed (dotted)). The (x) indicates the ACDM point (Qm = 
0.27, wo = — l,u)i = 0). The plain line separates accelerating (go < 
0) from decelerating (qo > 0) models. 

Simulation and interpretation : 

We have simulated, as in our previous paperQ, SNIa 
data corresponding to the same statistical power as the 
data sample, where we vary the fiducial values to study 
the effects of the priors (on Qm or/and wi). 

We start with some illustrations of the bias introduced 
by the J1m prior when it is different from the fiducial 
value. We consider two fiducial models which are com- 
patible with the data, when no prior is applied: one in 
acceleration with Q M =0.5, Wq =-2.2, wf =1.6 and one 
in deceleration with fl M =0.5, u>^=-0.6, wf =-10. We 
apply the 4-fit to the two models with the two strong 
priors: Q M = 0.27 ±0.04 and Q M = 0.5 ±0.04. FigureEl 
shows how the prior affects the conclusions: 
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FIG. 2: Fisher contours in the (u>o, wi) plane at 68.3%CL for the 
two fiducial models: (in acceleration) =0.5, Wq=-2.2, 10^=1.6 
and (in deceleration) f2^ =0.5, Wq =-0.6, w^=-l0. The upper 
indice F is added to avoid confusion between Fiducial values and 
fitted values. The plain big and small ellipses correspond to the 
first model analysed with the strong prior Qm = 0.5 ± 0.04 (big 
ellipse) or with fi M = 0.27 ± 0.04 as in (small ellipse). The 
dashed big and small ellipses correspond to the second model. 



• When the correct prior on VLm is applied, the central 
values are not biased but the errors are very large. 

• When the wrong prior VLm — 0.27 ±0.04 is applied, the 
fitted values are wrong but in agreement with the concor- 
dance model. The statistical errors are very small. In all 
cases, x 2 is good and does not indicate that something 
is wrong. 

• With the data, it is not possible to distinguish between 
these two models, but the prior value can lead to wrong 
conclusions both on values and errors of the fitted 
parameters. 

We have then performed a complete fit analysis on the 
simulated data and scanned a large plane of fiducial val- s 
ues (wQ,wf) with 3-fits and 4-fits, assuming a flat uni- 
verse and using two fiducial values for ; 0.27 or 0.5. 
We always use in the fitting procedures, the strong prior 
Qm = 0.27 ± 0.04. The case Om = f!^ is equivalent 
to a Fisher analysis and only the errors are studied. In 
the case Qjv^ ^ biases are introduced in the fitted 
values. 

Figure |3| shows the fitted wo and W\ iso-lines for the 
4-fits in the biased case. The iso-lines are straight lines 
(not shown on the figure) when f2w = 0.27 (unbiased 
correct prior), but are biased when f2^ = 0.5. This is 
due to the strong correlations between wq and fijvf, and 
between wq and wi . 

In this configuration, we observe that, for the 4-fit, 
when — 5 < Wq < (a relatively wide range), the fit- 
ted values for wq are in a narrow range centred on -1 : 
— 1.8 < wq < 0. For w\, the situation is even worse since 
with fiducial values —8 < wf < 8, we get essentially pos- 
itive values for the fitted w±. The actual shapes of the 
distortions between the fiducial and the fitted values are 



FIG. 3: Fitted wo (left) and w± (right) iso-lines with 4-fits for a 
fiducial with Qj^ =0.5 and a prior Clj^ = 0.27 ± 0.04, in the plane 
(tug ,w[ ). 

readable on Fig. |3| 

A similar analysis performed with the 3-fit shows that 
the situation is even worse : one gets —1.5 < w < 
whatever the value of Wq . As w\ is forced to and 
£Im to 0.27, wq is closer to -1 which corresponds to the 
preferred solution for the fit. 

One can illustrate further this very problematic point, 
by defining "confusion contours" , namely some contours 
which identify the models in the fiducial parameter space 
(e.g. (wq,Wi)) that could be confused with another 
model. For instance, the contours of Figure 0] give the 
models in the plane (wQ,wf) with f2f f = 0.5 that can 
be confused (at 1 and 2a) with the concordance model 
if the (wrong) strong prior is applied. The two models 
used for the illustrative Fig. [21 are taken from extreme 
positions in this confusion contour of Fig. 0] 

For the 3-fit, the confusion contours with the concor- 
dance model are very large and include all models having 
roughly wf < (~5wq — 10). The situation here is par- 
ticularly bad since the fitting procedure is making two 
strong assumptions = and £!m = 0.27±0.04) which 
are not verified by the fiducial cosmology (two biases). 
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FIG. 4: Confusion contours in the fiducial plane (wQ,wf ) which 
identify the models that would be confused with the concordance 
model at 1 and 2<r for the 3-fit(left) and the 4-fit(right) procedures. 
The strong prior is Qm = 0.27 ±0.04 whereas the fiducial model is 
= 0.5. The vertical line separates accelerating from decelerat- 
ing models. 



The next step is to study the parameter errors. We 
look at the correlation of the errors using fiducial models 
where £1^ =0.27 or 0.5. We determine the Wq and W\ 
errors, scanning the full plane (wq ,wf ) using 4-fits. 



Some regions of the parameter space (see Figure HJ) are 
favoured and always produce small errors. This is due to 
the correlation between wo and w\. The error depends 
strongly on the fitted wq and w\ values but not strongly 
on the VLm value: a different value of Q.m affects the scale 
of the errors but not the shape of the plots. We find a 
linear scaling of the error when we change Qm from 0.27 
to 0.5 (i.e. a factor 2). 
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FIG. 5: wq and u>i errors for a 4-fit with the correct strong prior, 
for a fiducial with fif^ =0.27. 

Combining this with the previous paragraph leads to 
an interesting point, i.e. the favoured fitted values of 
the fit (wq > —1.8 and wi > 0), which were shown to 
be mainly driven by the prior value, correspond also to 
the region of the plane where the parameters errors are 
always small. 

We conclude that the applied fitting procedure with 
this strong prior can bias the conclusions by constraining 
the (uio,wi) solution near the (-1,0) solution, where the 
statistical error is always very small. In particular, Riess 
et al.0 found a gain factor of order 8 on the accuracy 
of the neasurement compared to previous analysis. This 
is mainly due to the prior and not to the inclusion 
of the high z HST events. The present observational 
constraints on Qm are thus an important issue. 

Revisited conclusions on existing data : 

Most reviewers of cosmological parameters favour a 
value close to the strong prior choice made by 0- This 
result is based on WMAP[8j data combined with 2dF 
dataQ or more recently with SDSS data 10], and cor- 
responds to P(E3) &m = 0.27(0.3) ± 0.04 with h = 
0.71(0.70)1qq3. However, these results are based on sev- 
eral prior assumptions in order to lift the degeneracies 
among the various cosmological parameters (e.g. CIm, 
h, erg, w ...). For instance, Spergel et al.|j| mention that 
a solution with Q,m = 0.47, w = —0.5 and h — 0.57 in the 
CMB is degenerate with the ACDM model. This kind of 
solution is excluded for three reasons: the Hubble Con- 
stant value is 2ct lower than the HST Key Project value 
and the model is a poor fit to the 2dF and SN data. 

However, i) in spite of the precise HST result, the 
Hubble constant value is still controversial ii) 
We have shown in the previous section that the SN data 
analysis can only conclude that Am < 0.6 (see Fig. [Q. 
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FIG. 6: The shaded region is excluded at 95%CL when no prior 
on flu is applied. The ellipse corresponds to the strong prior 
constraints as in Q. The (x) is ACDM. The dotted line separates 
accelerating from decelerating models. 

Hi) The 2dF and the SDSS Collaborations [1 [H have 
extracted from an analysis of the power spectrum 

of galaxy redshift surveys. The degeneracy between D,Mh 
and the baryon fraction is lifted thanks mainly to the pre- 
cise determination of the baryon fraction by CMB data 
(see Fig.38 of 0). Should that prior change, the pre- 
ferred values from SDSS would indicate a higher value of 

In addition, a large variety of observations give con- 
straints on Hm, which is found to vary from 0.16 ± 0.05 
[3 up to a value above 0.85 [l5|. 

Conversi et al.|l2| provide an interesting critical anal- 
ysis on the present constraints on cosmological parame- 
ters, especially on Cljsti h, and w. Through the study of 
the degeneracies, they show that the result Qm — 0.27 ± 
0.04 is obtained under the assumption of the ACDM 
model, and provide specific examples with smaller h 
(h < 0.65) and higher Qm (Om > 0.35) which are in 
perfect agreement with the most recent CMB and galaxy 
redshift surveys. 

In conclusion, we follow the point of view of Bridle 
et al.0], who argue that it may be "that the real un- 
certainty is much greater" than the 0.04 error obtained 
from the combination of CMB and large scale structure 
data. 

Returning to SN data analysis, we suggest, for the time 
being, to reevaluate the conclusions by relaxing the £Im 
central value. Figure shows the 95%CL constraints 
in the (w a ,wi) plane obtained from the gold sample @ 
with no prior assumption on Qm ■ Taking an uncertainty 
of 0.2, which is the intrinsic sensitivity of SN results (see 
Table I, line 6), does not change the conclusions: 

• Large positive variations in time of the equation of state 
are excluded (at 95%CL) since the dark energy density 
blows up as e 3wiz [iH ]. 

• The quintessence models which have in general (wq > 
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— 1, 1 > loi > 0) are seriously constrained. For in- 
stance, the SUGRA model 0] characterized by w ~ 
—0.8 and w\ « 0.30 is close to the border of the 95%CL 
contour (precisely, one gets Ax 2 =3.5 corresponding to an 
exclusion at 80%CL). 

• The quadrant (wg > — 1, W\ < 0) corresponding to k- 
essence models 19] or some Big Crunch models pi l2fjj. 
is not the "least favoured" , contrary to the conclusions 
drawn with the strong prior Q. We find that if wq goes 
towards 0, then wi should be more and more negative. 

• If Wo < — 1) the constraints on wi are weak (except 
for large positive values). This region of the plane corre- 
sponds to phantom models [2lJ which have unusual prop- 
erties and may have very different consequences for the 
fate of the Universe (e.g. models with w± > wq + 1 will 
end in a Big Rip Models with very exotic w(z) may 
come from modified gravity!^ ■ The class of models with 
W\ < is roughly excluded at 95%CL, if the strong prior 
Q.M = 0.27 ± 0.04 is used Q, but is perfectly allowed for 
higher Qm values (or larger prior errors). 

• As can be seen on Fig. 6 (and also on Figures 1, 4 
and the decelerating model used to draw Fig. 2), our 
analysis without assumptions on and w\, allows 
decelerating models with specific properties : low wq, 
£Im ~ 0.5 and w± << 0. One can wonder if this result 
is not in tension with the geometrical test performed 
in where the only assumption is to use a linear 
functional form for q(z) (i.e. qo + q\z). It can be shown 



that a varying equation of state implies a non-linear 
q(z) 1 in particular, the linear approximation breaks 
down if w[ < —1.5. More details on this more sub- 
tle analysis will be presented in a forthcoming paper |23| . 

To go further, a coherent combined analysis of all data 
is mandatory, with a proper treatment of correlations 
and no prior assumptions. Some recent papers go in that 
direction [3 Hi El • 

In addition, as soon as the statistical errors will 
become smaller, systematic questions cannot be ne- 
glected and should be controlled at the same level of 
precision. This is the challenge for the next generation of 
experiments. A promising approach is to combine SNIa 
with weak lensing, as prop osed by the future dedicated 
SNAP/JDEM mission HJ. 
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